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Abstract The genus Rhabdophis is a group of widely 
distributed snakes with more than 20 species. Recent field 
surveys uncovered a species in southwestern China, which 
has long been considered as R. pentasupralabialis. Combined 
molecular and morphological analyses revealed it as а new 
species Rhabdophis chiwen sp. nov. Based on 12 specimens, 
this new species is distinguished by the following 
characters: 1) dorsal body saddlebrown, dorsal scales 
typically with black margins forming spots and stripes, the 
margin of the outer row forming two faint dorsolateral 
black cross-bars alongside body; 2) ventral scales 151-159, 
the outer margin of ventral scales and several lateral 
rows of dorsal scales forming ventrolateral longitudinal 
brownish-red coloration, with faint black spots in the 
middle of ventral scales; 3) a black oblique stripe present 
below eyes, often with a black spot between the 2nd and 
3rd supralabial and a black stripe on the 5th suprala bial; 4) 
eyes dark khaki, pupils black; 5) infralabials usually 7, the 
first four in contact with anterior chin-shields; 6) temporal 
scales 1+1; 7) dorsal scales in 15 rows, feebly keeled except 
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the outer 1-2 rows; 8) anal scale divided; subcaudals 45-59; 
9) preocular 1 and postoculars 3 (occasionally 2); 10) body 
medium-sized (snout-vent length: adult males 404-431 mm, 
adult females 409-476 mm); 11) tail moderate (tail length/ 
total length in adult males 0.205–0.238, in adult females 
0.172-0.193). With the discovery of this new species, the 
total number of species in genus Rhabdophis is 28 with 12th 
species known to occur in China. 


Keywords Rhabdophis chiwen sp. nov., morphology, 
Natricinae, phylogenetics, taxonomy 


1. Introduction 


The genus Rhabdophis is a group of widely distributed snakes 
with a unique defensive system that relies on nuchal glands 
(Mori et al., 2012). In 1960, Malnate defined this genus based 
on morphological characters, ie. hemipenes with sulcus 
spermaticus divided; last two maxillary teeth strongly enlarged, 
recurved and usually preceded by a diastema; internasals broad 
anteriorly, nostrils lateral; apical pits present or absent; vertebral 
glands present in several species (Malnate, 1960). 

Currently, 27 species are recognized and are distributed 
in eastern and southern Asia (Zhu et al., 2014; Takeuchi et 
al, 2018). Two new species were discovered іп the past five 
years, Rhabdophis akraios (Doria et al., 2013) and Rhabdophis 
guangdongensis (Zhu et al., 2014). 
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Eleven species of Rhabdophis are now known to occur 
in China (Zhu et al., 2014), which are R. adleri ZHAO, 1997, 
R. guangdongensis Zhu, Wang, Takeuchi and Zhao, 2014, R. 
himalayanus Günther, 1864, R. lateralis Berthold, 1859, R. leonardi 
Wall, 1923, R. nigrocinctus Blyth, 1855, R. nuchalis Boulenger, 1891, 
R. pentasupralabialis Jiang and Zhao, 1983, R. subminiatus Schlegel, 
1837, R. formosanus Maki, 1931 and R. swinhonis Günther, 1868. 
During several field work in the southwestern part of Sichuan 
Province, China, we collected a number of samples which were 
originally referred to as R. pentasupralabialis. Nevertheless, those 
snakes can be distinguished from R. pentasupralabialis based on 
morphological characters. Moreover, the phylogenetic analysis 
based on the mitochondrial cytochrome b (cyt b) and nuclear 
oocyte maturation factor (cmos) gene sequences revealed that 
those samples are different from R. pentasupralabialis and other 
congeners in China. Therefore, we propose that those specimens 


be a new species in the genus Rhabdophis. 
2. Materials and Methods 


2.1. Specimens sampling Specimens were collected from 
the Xingou Village (XG, 10 specimens) in Tianquan County 
and Jiguan Mountain (JG, 2 specimens) in Chongzhou City, 
Sichuan Province, China (Table 1). Fresh liver or muscle tissues 
were taken and immediately preserved in 95% ethanol and 
stored at –20 °C. The specimens were dehydrated in absolute 
ethanol and then later transferred to and stored in 70% ethanol. 
In addition, another 10 samples of R. leonardi (4 specimens), R. 
nuchalis (4 specimens) and R. pentasupralabialis (2 specimens) 
that represent 3 sympatric species with relatively close niche in 


Sichuan Province, were used for molecular analyses. 


2.2. Molecular phylogenetic analysis Genomic DNA was 
extracted from the collected tissue samples in Chengdu Institute 


Table 1 Morphological characters of Rhabdophis chiwen sp. nov. 
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of Biology (CIB). PCR amplifications were performed in 25 
pl reactions (12.5 pl I-5™ 2xHigh-Fidelity Master Міх, 10 pl 
ddH,O, 1 pl F-primers, 1 ul R-primers, 0.5 р DNA template) 
using the following cycling conditions: initial denaturation 
for 2 min at 95°C, followed by 35 cycles: denaturation at 94 
°C for 40 s, annealing at different temperatures (48.5 °C for cyt 
b or 56 °C for cmos) for 25 s, elongation at 72 °C for 15 s, and 
then finalized with elongation step of 2 min at 72 °C, with the 
PTC-100 thermal cycler (BioRad, USA). Purified PCR products 
were sequenced using the same PCR primers. Sequencing was 
completed by Beijing Qingke New Industry Biotechnology 
Co., Ltd. Sequences for comparison of available species were 
downloaded from GenBank. Amphiesma stolatum and Natrix 
natrix were used as outgroups. 

The cyt b and cmos sequences were combined into one dataset 
to build gene trees. All sequences were aligned with other 
retrieved sequences in the same gene loci by using software 
MEGA 7 (Kumar et al., 2016). Raw trace files were edited in 
Geneious 7 (Biomatters Limited, New Zealand). Partition 
finder 2.1.1 under BIC identified the optimal models of sequence 
evolution for each partition (Lanfear et al., 2012). Phylogenetic 
relationships derived from the combined gene fragments were 
performed based on Bayesian Inference (BI) by using MrBayes 
3.2 (Ronquist et al., 2012). We ran our analysis for 20 million 
generations with the chains, sampling every 1000 generations. 
The average standard deviation of split frequencies (ASDSF < 
0.01). The first 1000 trees (of 20,000) were discarded as the burn- 
in. A 50 majority-rule with compatible groups consensus was 
taken from the remaining trees and posterior probabilities 
(pp) of 0.95 or above were considered significant. We also 
performed maximum-likelihood (ML) analysis by using the 
program RaxML v8 (Stamatakis, 2014) and IQ tree (Nguyen et 
al, 2015). Reliability of the ML tree was assessed by calculating 
bootstrap probability (BP) with 1000 replications. Additionally, 


No. SVL Tal SPO РКО PTO LR TEM SL IL VEN DSR 5С Sex Locality 
Holotype CIB116092 422 118 1 1 LARES il ЛЕВИ 5 L:7R:6 152 15-15-15 56 Male JG 
Paratype CIB116093 431 131 1 3 1 141 5 7 155 5-15-15 57 Male JG 
Paratype CIB116094 411 Ш 1 1 3 i iil 5 7 154 15-15-15 55 Маје XG 
Paratype CIB116095 461 96 1 L:3R:2 1 1+1 5 L:6R:7 157 5-15-15 45 Female XG 
CIBDL1807132 427 128 1 1 3 il ile 5 L:8R:7 154 15-15-15 58 Male XG 
CIBDL1807133 409 128 1 3 1 141 5 7 157 5-15-15 58 Male XG 
CIBDL1807134 411 Oil 1 IBS ING il Gl 5 706201088159. 15-15-15 49 Female XG 
CIBDL1807135 476 107 1 3 1 141 5 7 156 5-15-15 46 Female XG 
CIBDL1807137 409 So i 1 3 il ier 5) 17 К:6 159 15-15-15 52 Female XG 
CIBDL1807138 427 122 L:i2R:1 ЛЕЗ 1 141 5 7 152 5-15-15 59 Маје XG 
CIBDL1807139 404 104 1 1 9 Т 5 7 151 15-15-15 55 Male XG 
CIBDL18071311 422 91 1 3 1 141 5 7 158 5-15-15 47 Female XG 


Note: See Material and Methods section for the abbreviations. 
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the uncorrected cyt b and cmos p-distance matrix was compared 
using MEGA 7.0 (Kumar et al, 2016). 


2.3. Morphological analysis Twelve adult snakes (7 males and 
5 females) were examined. Information on some morphological 
characters of known species were obtained from literature 
(Boie, 1827; Boulenger, 1896, 1900, 1906; Bourret, 1935; Cantor, 
1839; Das, 2010; David and Vogel, 2010; de Lang and Vogel, 
2006; Doria et al, 2013; Duméril et ај, 1854; Günther, 1858, 1864; 
Leviton, 1970; Smith, 1943; Stuebing and Lian, 2002; Tweedie, 
1953; Zhao, 1997; Zhao et al., 1998; Zhao and Adler, 1993; Zhao 
and Jiang, 1981) and are shown in Table 2. Color description 
was done according to wiki color-coding, 

Terminology for morphological measurements and 
descriptions is as follows: body and tail length were measured 
with a tape ruler to the nearest 1 mm; total length (TL), from 
the tip of snout to the tip of tail; snout-vent length (SVL), 
from the tip of snout to posterior margin of cloaca; tail length 
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(TaL), from posterior margin of cloaca to the tip of tail. Other 
measurements were conducted with a digital caliper to the 
nearest 0.1 mm: head length (HL), from the snout tip to the 
posterior margin of the mandible; head width (HW) was 
measured at the widest part of the head on posterior side; the 
eye horizontal diameter (ED), the greatest diameter of the orbit. 
Biometric measures were performed exclusively on the right 
eye. Ratio of tail length to total length (TaL/TL) was recorded. 
The dorsal scale rows (DSR) were counted at one head length 
behind head, at mid-body, and at one head length before vent; 
the number of ventral scales (VEN) was counted according to 
the method proposed by Dowling (Dowling, 1951), half ventral 
was counted as one. For subcaudals (SC), first scale under the 
tail meeting its opposite was regarded as the first subcaudal 
scale, and the unpaired terminal scute was not included in the 
number of subcaudals. Paired scales on head were counted 
on both sides of the head and presented in left/right order. 


Table 2 Significant characters of R. chiwen sp. nov. and the other known 27 species of the genus Rhabdophis. 


Species DSR Outer dorsal scale row VEN SC PRO PTO MTA 

R. chiwen sp. nov. 15-15-15 smooth 151-159 45-59 1-2 2-3 no data 

R. pentasupralabialis 5-15-15 smooth 135-162 43-64 2-3 18 

R. adleri 19-19-17 feebly keeled 150-164 73-88 1-2 3-4 27(25+2) 

R. akraios 9-19-17 keeled 167-184 44-66 3 19(17+2) 

R. angeli 16(15)-15-15 smooth 117-126 39-46 1 3 22-23 

R. auriculata 7-17-15 strongly keeled 157-162 76-91 -2 3 27-32 

R. barbouri 19 at midbody strongly keeled 148-166 96-101 2 3 no data 

R. callichroma 9 at midbody keeled 152-159 79-86 -2 3 27-35 

R. callistus 21 at midbody strongly keeled 156 76 1 4 no data 

R. ceylonensis 9 at midbody strongly keeled 137-143 40-54 3 24–26(22–24+2) 
R. chrysargoides 21 at midbody strongly keeled 154-161 64-79 1 3 24(22+2) 

R. chrysargos 9 at midbody strongly keeled 143-175 60-93 -2 3 27-35 

R. conspicillatus 19 at midbody keeled or smooth 138-147 40-53 1 3 no data 

R. flaviceps 9 at midbody keeled 120-138 49-60 3-4 22-23 

К. guangdongensis 15-15-15 smooth 126 39 1 2 20 

R. himalayanus 9-19-17 feebly keeled 165-171 82-88 2 3 26(24+2) 

К. leonardi 18(17)-17-15 smooth 149-159 43-62 1 2=3 19(17+2) 

К. lineatus 9 at midbody strongly keeled 132-142 66-71 2 3 18(16+2) 

R. murudensis 19-19-15(17) feebly keeled 176-185 63-97 1 3 23(21+2) 

R. nigrocinctus 9-19-17 keeled or smooth 150-170 80-97 1 3-4 28(26+2) 

К. nuchalis 15-15-15 smooth 144-169 35-65 1 3 18-22 

R. plumbicolor 23-27 at midbody strongly keeled 144-160 27-45 2 3-4 no data 

R. rhodomelas 19 at midbody strongly keeled 126-136 41-56 1-2 3-4 no data 

R. rudis 23-23-19 strongly keeled 123-155 44-60 3 3-4 13-20(11-18+2) 
R. spilogaster 19 at midbody feebly keeled 148-156 75-92 1-2 3-4 15 

R. subminiatus 19-19-17 keeled or smooth 144-184 56-97 1 3-4 23-26(21-24+2) 
R. swinhonis 15-15-15 weakly keeled or smooth 124-165 44-74 1 2-3 19-23(17-21+2-3) 
R. tigrinus 19-19-17(15) keeled or smooth 144-188 38-74 1-2 2-4 22-23(20-21+2) 


Notes: See Material and methods section for the data sources. 


A : total number of maxillary teeth (the number of maxillary teeth before the diastema + the number of maxillary teeth behind the diastema). 
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Supralabials (SL) were considered being those shields that were 
behind the rostral and bordering the mouth gap; infralabials 
(IL) were considered being those shields that were behind the 
mental, completely below a supralabial and bordering the 
mouth gap. Number of supraocular (SPO), preoculars (PRO) 
and postoculars (PTO) were counted as shields above, at the 
anterior and posterior margin of the orbit. Loreal scales (LR) 
were counted as scales between the nasal scale and preocular. 
Number of temporal scales (TEM) were counted as scales behind 
the postoculars and between the parietal scales and supralabials. 
The sex was determined by inspection of the existence of 
hemipenes. 

This new species is most similar to and has long been 
considered as R. pentasupralabialis. To illuminate this problem, 
22 specimens (11 males and 11 females) collected from the 
type locality (Jiulong County, JL) of R. pentasupralabialis were 
examined and measured. Morphological characters of specimens 
were analyzed using Principal Component Analysis (PCA) 
(Wiister et al., 1992). Male specimens were from Jiulong County, 
Jiguan Mountain and Xingou Village, using the characters 
VEN, SC, left IL, right IL, left PRO, left PTO and right PTO. 
Female specimens were from Jiulong County and Xingou 
Village, using the characters VEN, SC, left IL, right IL, left PTO 
and right PTO. All statistical analysis were conducted with 
SPSS 25.0 (IBM Inc, Armonk, US.A,). 


3. Results 


31. Molecular phylogenetic analysis We obtained alignments 
of the mitochondrial gene cyt b (1074 bp) and nuclear gene cmos 
(535 bp). Sequence data were uploaded to GenBank, available 
accession numbers showed in supplementary Table 51. The best 
evolution models of each partition combination are shown in 
Table 3. The result show that cyt b and cmos genes consistent 
phylogenetic trees were achieved by using Bayesian Inference 
(BI) and Maximum likelihood RaxML (Figure 1). 

The topological structures of combined cyt b and cmos 
sequences phylogenetic tree are identical with the earlier study 
(Takeuchi et al., 2018). Our results also strongly supported the 
monophyly of Rhabdophis. The individuals from the Jiguan 
Mountain and Xingou Village were placed in the genus 
Rhabdophis with a high support (Figure 1). The Results also 
showed that they share a common ancestor with R. nuchalis, R. 
leonardi, and R. pentasupralabialis. 

The specimen (HK V36838) from Hongya County, Sichuan, 
China was also nested within the clade with individuals from 
Jiguan Mountain and Xingou Village. 

The uncorrected p-distances between species are shown in 
Table 4. The pairwise distances between species in the genus 
were wide, ranging from 0.004 to 0.198. The largest disparity 


| Vol. 11 


Table 3 The best evolution models of each partition combination. 


Subset Best Model Site (bp) Subset Partition 

1 HKY+I+G 358 cytb 1st 

2 HKY+I+G 554 cmos 3rd, cyt b 2nd 
3 GTR+G 358 cyt b 3rd 

4 K80+I 392 cmos 1st, cmos 2nd 


occurred between R. subminiatus from Thailand and R. 
chrysargos from Malaysia (0.198). In comparison, the specimen of 
Rhabdophis from Jiguan Mountain (CIB116092) differed from R. 
pentasupralabialis from Jiulong County, Sichuan, China (GP1065, 
type locality) by 0.081, and was different from R. guangdongensis 
(SYSr000018) by 0.066 substitutions per site to the other species. 
Moreover, the specimen (HK V36838) from Hongya County, 
Sichuan, China was identified as R. nuchalis before (Alfaro and 
Arnold, 2001), but its pairwise distance with R. nuchalis (HT0854 
and CIBDL1807208) were rather large, by 0.063 and 0.058 
respectively. 


3.2. Taxonomic conclusions The individuals from the Jiguan 
Mountain and Xingou Village were placed as monophyletic 
clade with strong supports (100% in PP and BP). In addition, 
the p-distance between the individuals of the clade and other 
species indicated very large values (more than 0.058). Combined 
with further evidence from morphology mentioned below, the 
individuals from southwestern Sichuan Province represent an 
undescribed species of the genus Rhabdophis. Here, we described 
it as a new species. 


Taxonomy 

Rhabdophis chiwen sp. nov. Chen, Ding, Chen and Piao, 2019 
(Figures 2, 3, 4, 6); 

Rhabdophis nuchalis pentasupralabialis: Jiang and Zhao, 1983 pp. 
59-62; 

Rhabdophis pentasupralabialis: Zhao, 1998, pp. 271-274; Zhao, 2006, 
pp.268-269; Takeuchi et al, 2018, Figure 2 and 3, p. 10226 (part). 


Diagnosis. 1) nuchal groove present, with enlarged and paired 
scales on each side; 2) Dorsal body saddlebrown, DSR in 15 
rows throughout, feebly keeled, the outer 1-2 rows smooth; 3) 
dorsal scales typically with black margins forming some spots 
and stripes, the margin of the outer row forming two faint 
dorsolateral black cross-bars alongside body; 4) a black oblique 
stripe below the eye, often with a black spot between the 2nd 
and 3rd SL and а black stripe (or separated as black spots) on the 
5th SL and the Ist TEM; 5) Eyes dark khaki, pupils black; 6) SL 
5, the 3rd and 4th touching the eye: 7) usually 7 ILs (occasionally 
6 or 8), the first four contact with anterior chin-shields; 8) 1 LR, 
and TEM 1+ђ 9) VEN 151-159, the outer margin of ventral scales 
and several lateral rows of dorsal scales forming ventrolateral 
longitudinal brownish-red coloration, with faint black spots in 
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Figure 1 Phylogenetic tree derived from the combined gene fragment (cmos + cyt b). Bayesian posterior probabilities (PP) and bootstrap 


probabilities (BP) from maximum likelihood are shown at the nodes. 
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the middle of each ventral scale and between scales which line a 
black stripe; 10) anal divided; SC 45-59; 11) PRO 1 (rarely 2) and 
PTO usually 3 (occasionally 2); 12) medium-sized body (SVL of 
adult males: 404-431 mm and adult females: 409-476 mm); 13) 
tail moderate and longer in males than in females (adult males: 
104-131 mm, adult females: 91-107 mm). 


Etymology. The species name of the new species “chiwen” is in 
reference to the ninth son of Loong in ancient Chinese myth 
who likes eating fire, and indicates the firefly-eating habit 
of this new species (Yoshida et al., 2020). Its common name is 
suggested as “Chiwen Keelback” in English and “ ИЕ” 
in Chinese. 


Holotype. CIB116092, an adult male, collected from Jiguan 
Mountain, Chongzhou City, Sichuan Province (30°4557.64" N, 
103°143.52' E, and 1846 m as.) collected on Ist July 2019 by Li 
Ding (Figure 2 and 3). 


Paratypes. СІВ116093, an adult male, collected from the same 
locality as holotype on Ist July 2019 by Li Ding. CIB116094- 
95, two specimens, an adult male and an adult female, collected 
from Xingou Village, Tianquan County, Sichuan Province 
(29°55:42.39" N, 102°239.19" E, and 1461 m asl) collected on 22nd 
July 2018 by Li Ding and Zening Chen. 


3.3. Description of the holotype Adult male with TL 540 
mm (SVL 422 mm and TaL 118 mm). TaL/TL 21.9%. Body 
elongated and cylindrical. Nuchal groove present; seven scales 
on each side of groove more or less distinctly enlarged and 
paired. DSR in 15 rows throughout, feebly keeled, the outer 2 
rows smooth. VEN 152; anal divided; SC 56, paired, final spine 
present at tip of tail. 

Head oval and distinctly wider than neck, HL 13.67 mm, 
2.53% of TL; HW 9.10 mm, HL/HW 1.50. Eyes are relatively 
large, with rounded pupil, EW 2.95 mm, EW/HL 21.58%. Rostral 
semi-circular when viewed from the front, wider than long, 
barely visible from above. Nostrils large and open laterally, 
nasal completely divided, in contact with rostral. Internasals 
slightly shorter than prefrontals. Prefrontals wide, bending 
to the loreal region. Frontal slightly pentagonal, longer than 
wide, equal to the distance from the rostral to frontal. LR 1, 
longitudinally longer. SPO 1; PRO 1; PTO left 2 and right 3; 
subocular absent. SL 5, the 3rd and 4th touching the eye, with 
a black oblique stripe below the eye, a black spot between the 
2nd and 3rd SL and a black stripe on the 5th SL and the Ist 
TEM. Mental triangular, width approximately twice of length. 
IL left 7 and right 6, the first one in contact with each other 
behind the mental and the first four contact with anterior chin- 
shields, the 4th to 6th infralabials contact with posterior chin- 
shields. At the left side, the 7th infralabial slightly longer than 
the 6th infralabial, but much narrower, width about 1/3 of its 
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Figure 2 Dorsal (A) and ventral (B) view of the holotype of Rhab- 
dophis chiwen sp. nov. (CIB116092). Photo by Shengchao Shi. 


height. TEM 2 each side, anterior temporal single and posterior 
temporal single. 

Hemipenis bilobed, extends to the 11th subcaudal, forked 
at the level of the 9th subcaudal. When wholly everted, the 
hemipenes have small and hard spines resembling hooks 
covering the entire organ with the highest density of spines on 
the tip. The length of each lobe is about 1/4 of the hemipenis. 
Basal big spine absent, skin shallow cupped. The centripetal 
sulcus bifurcates at the fork and extends to the tip of per lobe. 
Lips evident and smooth (Figure 4). 


Coloration of specimens in life. Dorsal body saddlebrown 
with rusty red. Dorsal scales are typically with black margins 
forming some scattered spots and stripes (for holotype, those 
dispersed black spots forming a black stripe around the neck, 
see Figure 3 A, C, E), the margin of the outer row forming two 
faint dorsolateral black cross-bars alongside body (Figure 2 A; 
Figure 5 A, C). Coloration of ventral is seashell and the outer 
margin of ventral scales and several lateral rows of dorsal scales 
forming ventrolateral longitudinal brownish-red coloration, 
with faint black spots in the middle of each ventral scale and 
between scales at the anterior part of body which line a black 


stripe and merge into black patches covering the whole scale 
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Figure 3 The holotype of Rhabdophis chiwen sp. nov. (CIB116092) in life. A-E: Different views of head. С and Е: showing the same 5 su- 
pralabials as R. pentasupralabialis. D: ventral view of head, displaying the representative 7 infralabials. F: Lateral view of trunk, demonstrat- 


ing several smooth outer rows. Photo by Shengchao Shi. 


with dispersed brownish-red spots at the posterior part of body 
(Figure 2 B). Hatchlings similar to adults except for darker 
dorsal coloration, a yellow stripe present on nape, separated by 
the nuchal groove (Figure 5 D); ventral scales typically shallow 
black. 


Coloration of specimens in preservative. Coloration pattern 
in preservative similar to alive, but the faint black spots in 
the middle of each ventral scale may become much dispersed, 
forming large black patches alongside the ventral surface after 
one year in preservative. 


3.4 Variation Other specimens generally resemble the 
holotype except the following characters (Table 1): IL rarely 8. 
PRO rarely 2 and PTO occasional 2. SVL ranges from 409-476 
mm in females and 404-431 mm in males. Female specimens 
tend to have a shorter tail and fewer SC than male specimens 
(for TaL 91-107 mm versus 104-131 mm, for SC 45-52 versus 
55-59). The ventrolateral longitudinal brownish-red coloration 
may be lighter (Figure 5 A, C) and faint black spots on the 
ventral scales could merge into large patches and cover the 


whole scale at the posterior part of body (Figure 5 B). The dorsal 
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Figure 4 The hemipenis of the holotype of Rhabdophis chiwen sp. nov. (CIB116092). Photo by Shengchao Shi. 
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Figure 5 Photos of Rhabdophis chiwen sp. nov.. A and B: Ventral view of an adult female displaying the coloration variations, photo by 
Yige Piao. C: An adult female in life. D: A yearling in life. For C and D, photo by Masaya Fukuda. 
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scales of some specimens (such as СІВІ16093) are feebly keeled 
even on the second outer row, with only the outer row smooth. 


3.5. Comparisons Comparative data of R. chiwen with 27 
known species of the genus Rhabdophis were obtained from 
literature (see Materials and methods) and are shown in Table 2. 

Currently, 11 species are known to occur in China, and 
Rhabdophis chiwen can be distinguished from these species 
by following characters (a detailed comparison between R. 
chiwen and R. pentasupralabialis is separately listed in a latter 
paragraph). R. chiwen differs from R. adleri and R. himalayanus 
by DSR 15 rows throughout vs. DSR 19-19-17 in R. adleri and 
R. himalayanus. R. chiwen differs from R. formosanus and R. 
lateralis by DSR 15 rows throughout vs. DSR 19-19-17(15) in 
R. formosanus and R. lateralis. The new species differs from 
R. guangdongensis by 5 SLs, VEN 151-159, SC 45-59, dorsal 
body saddlebrown vs. 6 SLs, VEN 126, SC 39, brownish-grey 
coloration in R. guangdongensis. The new species is different from 
R. leonardi by DSR 15 rows throughout vs. DSR 18(17)-17-15 in 
R. leonardi. The new species is different from R. nigrocinctus and 
R. subminiatus by DSR 15 rows throughout vs. DSR 19-19-17 in 
R. nigrocinctus and R. subminiatus. The new species differs from R. 
nuchalis by 5 SLs, 7(6,8) ILs, TEM 1+1, dorsal body saddlebrown 
vs. 6 SLs, 8(7) ILs, TEM 1+2, olive green coloration with black 
and magneta spots іп R. nuchalis (Figure 6). R. chiwen differs from 
R. swinhonis by 5 SLs, 7(6) ILs, dorsal body saddlebrown vs. 6 SLs, 
8(7) ILs, medium brown coloration with several rows of black 
spots in R. swinhonis. 

Rhabdophis chiwen, R. pentasupralabialis, R. guangdongensis, R. 
swinhonis, R. angeli and R. nuchalis have the minimal number of 
dorsal scale rows in this genus, equal to 15 rows. Nevertheless, R. 
chiwen can be readily distinguished from the other five species 
and all other species of Rhabdophis by SLs, DSR, VEN, SC, 
PRO, PTO or different coloration. Because R. guangdongensis, 
R. swinhonis and R. nuchalis are distributed in China, their 
comparisons with R. chiwen are listed in the previous paragraph. 
The new species differs from R. angeli by VEN 151-159, SC 45-59, 
dorsal body saddlebrown vs. VEN 117-126, SC 39-46, brownish 
coloration with a dorsolateral series of small reddish spots in R. 
angeli. 

For the remaining species within the genus Rhabdophis 
except R. pentasupralabialis, R. chiwen sp. nov. can be readily 
distinguished by its dorsal scales in 15 rows throughout vs. 
DSR 17-17-15 in R. auriculata; DSR 19 at midbody in R. barbouri, 
К. callichromus, R. chrysargos, R. conspicillatus, R. lineatus and R. 
spilogaster; DSR 21 at midbody in R. callistus and R. chrysargoides, 
DSR 19-19-15(17) in R. murudensis; DSR 19-19-17(15) in R. tigrinus. 

Specimens collected from Xingou Village and Jiguan 
Mountain differ from specimens of R. pentasupralabialis by: 1) 
larger size (mean TL 536.83 mm) vs. mean TL 483.64 mm in R. 
pentasupralabialis (Table $2); 2) saddlebrown coloration of dorsal 
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body vs. dark green or olive green coloration of dorsal body 
in R. pentasupralabialis; 3) the margin of the outer row forming 
two faint dorsolateral black cross-bars alongside body vs. absent 
in R. pentasupralabialis; 4) the outer margins of ventral scales 
and several lateral rows of dorsal scales forming ventrolateral 
longitudinal brownish-red coloration, with faint black spots in 
the middle of each ventral scale and between scales which line a 
black stripe vs. absent in R. pentasupralabialis (for the comparison 
of ventral coloration, see Figure 51). 

For male specimens, the first two principal components 
accounted for 55.900% of cumulative coefficients. Principal 
component 1 (РСІ) accounted for 31.464% and principal 
component 2 (PC2) for 24.436%. For female specimens, the first 
two principal components accounted for 64.567%. РСТ accounted 
for 39.106% and PC2 for 25.461%. Table 5 displays which 
characters are important in PC1 and PC2 for male and female 


specimens. 


Figure 7 shows the plots of the first two principal components 


A 


Figure 6 The lateral or ventral view of R. pentasupralabialis (A), 
holotype (CIB116092) of R. chiwen sp. nov. (B) and R. nuchalis (C), 
displaying the 6 ILs of R. pentasupralabialis, 7 ILs of R. chiwen sp. 
nov. and 8 ILs of R. nuchalis. 
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Table 5 Component matrix of Principal Component Analysis for 
male and female specimens. 


Component matrix 


Males Females 
Component Component 

1 2 1 2 
VEN -0.759 -0.012 0.697 0.543 
SC -0.351 0.452 -0.452 0.557 
left IL 0.647 0.038 0.502 -0.378 
right IL 0.781 –0.069 0.487 -0.756 

left PRO 0.395 0.821 - - 
left PTO -0.024 -0.860 0.785 0.139 
right PTO 0.563 -0.295 0.742 0.435 


Note: See Material and methods section for the abbreviations. 


for males (A) and females (B). The PCA for male specimens did 
not reveal much differentiation among the localities, while the 


female specimens could be easily distinguished from each other. 


3.6. Distribution, habitat and behavior Rhabdophis chiwen is 
currently known to be distributed in several parts of Sichuan 
Province, including Xingou Village of Tianquan County, 
Jiguan Mountain of Chongzhou City and Hongya County. 
During several field surveys in Xingou Village, individuals of 
this species were commonly encountered, including hatchings, 
juveniles and adults. At more than six different field sites with 
the GPS records, we speculate that Rhabdophis chiwen lives at 
the altitude range of 1100-2200m, typically near farmland and 
the source of water (Figure 8). With the stomach content by 
forced regurgitation and observation of feeding behaviors in 
the laboratory, it has been confirmed that this species primarily 
prey on earthworms and fireflies (Yoshida et al., 2020). The 
defensive behavior of this species is the typical body lift 
described by Mori (Mori et al., 2016). 


PCA for male specimens 
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4. Discussion 


The genus Rhabdophis is widely distributed among Eurasia. 
Species of this genus are known to occupy a wide range of 
microhabitat and different altitudes (Zhao, 2006). Sichuan 
Province lies in the southwestern part of China, which stretches 
across Palearctic Realm and Oriental Realm, and can mainly 
be divided into two parts: the western part is the West Sichuan 
High Plateau and the eastern part is the Sichuan Basin (Zhao, 
2002). The complex natural environment provides varied 
habitats for a large quantity of species (Wang et al., 2010). 
Currently, 5 species of Rhabdophis are reported to inhabit 
Sichuan Province, ie. R. pentasupralabialis, R. nuchalis, R. leonardi, 
R. subminiatus and R. lateralis. There are doubtful records of R. 
nigrocinctus, because this species typically inhabits low altitude 
habitat with high humidity rather than dry-hot valleys as the 
recorded locality in southwestern Sichuan (Li Ding, personal 
communication). R. chiwen was also found in Hongya County 
during several field surveys, where a specimen (HK V36838) was 
considered as R. nuchalis in a previous study (Alfaro and Arnold, 
2001). With our long-term field records since 2012, it is shown 
that the distribution range of R. chiwen is partly overlapped 
with that of R. pentasupralabialis and R. nuchalis. Together with 
К. leonardi, these four species all primarily feed on earthworms 
and form a earthworm-eating clade while most species in this 
genus prey on frogs and fish (Zhao, 2006; Yoshida et al., 2020). 
Intriguingly, all these earthworm-eating species can be found 
in Sichuan, and it remains to be investigated whether and how 
they show some niche differentia tions. 

A previous study has pointed out substantial genetic 
divergence within R. nigrocinctus, R. swinhonis, R. nuchalis, and 
especially R. subminiatus, and suggested the possibility for several 
undescribed species (Takeuchi et al., 2018). In the study of this 
genus, our field records clearly support this possibility (Li 
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Figure 7 The plots of the first two principal components for males (А) and females (В) of specimens. 
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Figure 8 Habitat of Rhabdophis chiwen in Xingou Village, Tian- 
quan County, Ya’an City, Sichuan Province. Photo by Yige Piao 
on June 2nd, 2019. 


Ding, personal communication). It is sure that the western and 
southern parts of Sichuan Province still provide shelters for 
certain cryptic species of Rhabdophis. 
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Figure S1 The ventral view of R. pentasupralabialis (CIB10668) (A) from type locality and R. chiwen sp. nov. (CIBDL1807 137) (B) from Xin- 
gou Village, displaying the black spots and patches on the ventral scales in R. chiwen sp. nov. vs. absent in R. pentasupralabialis. Both speci- 
mens were kept in preservative for more than one year. 


Table 51 Information of all sequence data used in the study. 


Sample ID Species Individual No. Country and Locality cyt b cmos 


2 Natrix natrix LSUMZ41506 Sheppey Island, Kent, United Kingdom AF471059 AF471121 


4 Rhabdophis lateralis_2 GP613 China GQ281785 JQ687444 


6 Rhabdophis tigrinus SCUM090004 Tianshui, Gansu, China KF765396 KF765391 


8 Rhabdophis subminiatus RDQ200905366 Baoshan, Yunnan, China KF765398 КЕ765393 


10 Rhabdophis subminiatus SICAU090029 Hong Kong, China KF800928 KF800918 


12 Rhabdophis guangdongensis SYSr000018 Aizhai, Renhuay, Guangdong, China KF800930 KF800920 


14 Rhabdophis leonardi RDQ200905367 Dali, Yunnan, China KF800932 KF800922 


16 Rhabdophis pentasupralabialis GP1065 Jiulong, Sichuan, China KF800934 KF800924 


18 Rhabdophis nigrocinctus CAS215280 Shan State, Myanmar KF800936 KF800926 


20 Rhabdophis himalayanus_1 HT0847 China LC325299 LC325746 


22 Rhabdophis nuchalis_4 HT0854 China LC325301 LC325748 


24 Rhabdophis swinhonis_1 HT0021 Taiwan, China LC325303 LC325749 


26 Rhabdophis tigrinus_1 HT0098 Japan LC325305 LC325751 


28 Rhabdophis nigrocinctus_ 1 HT0253 Thailand LC325307 1С325753 


30 Rhabdophis chrysargos HT0342 Malaysia LC325313 LC325759 


32 Rhabdophis subminiatus_2 HT0344 Thailand LC325315 LC325761 


34 Rhabdophis callichromus_1 HT0654 Vietnam LC325324 LC325770 


36 Rhabdophis subminiatus_4 HT0680 Vietnam LC325328 LC325774 


38 Rhabdophis pentasupralabialis_2 НТ0700 China LC325332 LC325778 


40 Rhabdophis swinhonis_2 НТ0717 Taiwan,China LC325334 LC325780 


42 Rhabdophis conspicilatus НТ0791 Malaysia LC325342 LC325788 


44 Rhabdophis nuchalis_3 HT0807 China LC325353 LC325799 


46 Rhabdophis adleri_ 1 HT0831 China LC325356 LC325802 


48 Rhabdophis nuchalis HKV36838 Hongya County, Sichuan, China AF402907 


(Continued Table 81) 


Sample ID Species Individual No. Country and Locality cyt b cmos 
50 Rhabdophis chiwen sp. nov. CIB DL1807203 Jiguan Mountain, Sichuan, China MN656320 MN656341 
51 Rhabdophis chiwen sp. nov. CIB DL1807204 Jiguan Mountain, Sichuan, China MN656321 MN656342 
52 Rhabdophis chiwen sp. nov. CIB DL1807206 Jiguan Mountain, Sichuan, China MN656322 MN656343 
53 Rhabdophis nuchalis CIB DLHK1 Chongqing, China MN656323 MN656344 
54 Rhabdophis chiwen sp. nov. CIB DL2014M1 Jiguan Mountain, Sichuan, China MN656324 MN656345 
55 Rhabdophis chiwen sp. nov. CIB DLM415 Jiguan Mountain, Sichuan, China MN656325 MN656346 
56 Rhabdophis nuchalis CIB DL2014M5 Longnan, Gansu, China MN656326 MN656347 
57 Rhabdophis leonardi CIB DL1806011 Panzhihua, Sichuan, China MN656327 MN656348 
58 Rhabdophis leonardi CIB DL1807031 Panzhihua, Sichuan, China MN656328 MN656349 
59 Rhabdophis pentasupralabialis CIB DL1807061 Panzhihua, Sichuan, China MN656329 MN656350 
60 Rhabdophis pentasupralabialis CIB DL1809267 Panzhihua, Sichuan, China MN656330 MN656351 
61 Rhabdophis nuchalis CIB DL1807207 Tangjiahe, Sichuan, China MN656333 MN656354 
62 Rhabdophis nuchalis CIB DL1807208 Tangjiahe, Sichuan, China MN656334 MN656355 
63 Rhabdophis nuchalis CIB DL1807209 Tangjiahe, Sichuan, China MN656335 MN656356 
64 Rhabdophis chiwen sp. nov. CIB116094 Xingou Village, Sichuan, China MN656336 MN656357 
65 Rhabdophis chiwen sp. nov. CIB DL1804151 Xingou Village, Sichuan, China MN656337 MN656358 
66 Rhabdophis chiwen sp. nov. CIB116095 Xingou Village, Sichuan, China MN656338 MN656359 
67 Rhabdophis leonardi CIB DL1809265 Yingjiang, Yunnan, China MN656318 MN656339 
68 Rhabdophis chiwen sp. nov. CIB116092 Jiguan Mountain, Sichuan, China MN656331 MN656352 
69 Rhabdophis chiwen sp. nov. CIB116093 Jiguan Mountain, Sichuan, China MN656332 MN656353 
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